Risk-sensitive Dissipativity of Linear Quantum Stochastic Systems 
under Lur'e Type Perturbations of Hamiltonians 



Igor G. Vladimirov, Ian R. Petersen 



Abstract 

This paper is concerned with a stochastic dissipativity theory using quadratic-exponential storage functions 
". for open quantum systems with canonically commuting dynamic variables governed by quantum stochastic 

differential equations. The system is linearly coupled to external boson fields and has a quadratic Hamiltonian 
which is perturbed by nonquadratic functions of linear combinations of system variables. Such perturbations are 
similar to those in the classical Lur'e systems and make the quantum dynamics nonlinear. We study their effect 
on the quantum expectation of the exponential of a positive definite quadratic form of the system variables. 
This allows conditions to be established for the risk-sensitive stochastic storage function of the quantum system 
to remain bounded, thus securing boundedness for the moments of system variables of arbitrary order. These 
results employ a noncommutative analogue of the Doleans-Dade exponential and a multivariate partial differential 
version of the Gronwall-Bellman lemma. 
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I. Introduction 

Qh- Quantum systems, originated in quantum physics and its applications, including, for example, the 
interaction of coherent light with matter at an atomic level studied in quantum optics [6], [30], can be 
described briefly as noncommutative stochastic systems and are treated by quantum probability tools. 
Their dynamic variables are represented in terms of an algebra of (generally, noncommuting) operators 
on a Hilbert space, with self-adjoint operators (usually referred to as observables) corresponding to real- 
valued physical quantities. The role of a probability measure is played, though in a noncommutative 
fashion, by a density operator [20] on the underlying Hilbert space which specifies the quantum state of 
the system. An efficient language to describe open quantum systems, interacting with the environment, is 
provided by quantum stochastic differential equations (QSDEs) [9], [23] (see also a review paper [7] and 
(*C) references therein) which govern the system variables in the Heisenberg picture of quantum dynamics. 
• The QSDEs are driven by a quantum noise from the surrounding memoryless heat bath, represented 
by a boson Fock space [23], and employ a system-bath interaction model in combination with the 
Hamiltonian which would specify the "internal" dynamics of the system in isolation from the external 
fields. This unified Markovian approach to open quantum systems and their interconnections is widely 
used in quantum control [1], [2], [4], [11], [12], [13] (see [34] for a more complete bibliography) which 
is aimed at synthesizing such systems and enhancing their performance (including stability, optimality, 
robustness) for practical applications, such as quantum metrology and quantum-optical communication, 
to mention a few [21]. An important class of open quantum systems is formed by those where the system 
variables satisfy canonical commutation relations (CCRs) (similarly to the quantum-mechanical position 
and momentum operators [20]), the system Hamiltonian is quadratic and the system-bath coupling 
operators are linear with respect to the variables. Such systems are governed by linear QSDEs and 
are dynamically equivalent to open quantum harmonic oscillators which are basic models in linear 
quantum stochastic control [4], [22], [25], [29]. Despite their theoretic convenience (including closed- 
form solutions and preservation of Gaussian nature of the system density operator [10], [25]), the linear 
quantum dynamics remain a simplified model, which may be distorted by an unmodelled energetics 
of the real system and its environment. The resulting perturbed behaviour can be "localized" in terms 
of bounds for the second moments of system variables, which play the role of storage functions to 
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quantify robust mean square stability [26] of the system in the quantum stochastic version [12] of the 
classical dissipativity theory [32]. In the present paper, we develop this line of research further by 
extending risk- sensitive storage functions [19] from classical stochastic systems to the quantum setting. 
The Hamiltonian of the quantum systems being considered has a nominal quadratic part which is 
perturbed by nonquadratic functions of linear combinations of the system variables. Such perturbations 
are similar to those in the classical Lur'e systems [16] and lead to nonlinear QSDEs which neither lend 
themselves to closed-form solution nor maintain Gaussian quantum statistics. We study the dynamics of 
the quantum expectation of the exponential of a positive definite quadratic form of the system variables. 
For this purpose, we use a noncommutative analogue of the Doleans-Dade exponential [3] based on 
a commutator approach to parameter differentiation of exponential operators [17], [18], [28], [31] to 
take advantage of the CCRs between the system variables. In combination with a multivariate partial 
differential version of the Gronwall-Bellman lemma, this allows conditions to be established for the 
perturbations of the Hamiltonian under which the risk- sensitive storage function of the quantum system 
remains bounded, thus securing boundedness of moments of the system variables of arbitrary order. 
The methods, developed in this paper, differ from those in [2], [11], [33] and can be of interest for the 
design of coherent (measurement-free) risk-sensitive quantum control schemes. 

II. Underlying quantum stochastic systems 

We consider an open quantum system with n dynamic variables X\(t), . . .,X n (t), which evolve in time 
t ^ and are assembled into a vector X[t) := (Xj(t))i^j^ n (vectors are organised as columns unless 
indicated otherwise). The system variables Xi(0),...,X n (0) at the initial moment of time are self- 
adjoint operators on a complex separable Hilbert space M 1 (such as the quantum-mechanical position 
and momentum operators) which are assumed to satisfy CCRs 

[X,X T ] := ([Xj,X k ]) 1<jMn =XX T -(XX T ) T = i&. (1) 

Here, [A, B] := AB — BA is the commutator of operators, and the transpose (-) T applies to matrices 
with operator- valued entries as if the latter were scalars. Also, i := y/— 1 is the imaginary unit, and 
:= (Qjk)i^j,k^n is a rea l antisymmetric matrix of order n (the space of such matrices is denoted by 
A w ). By the standard convention, a linear operator ^ on jtf? is lifted to its ampliation t, Cg)/jr on the tensor 
product space Jtf?®^, where /jr denotes the identity operator on a boson Fock space & [23]. The latter 
provides a domain for the action of an m-dimensional quantum Wiener process W(t) := (Wk(t))\^k^m 
adapted to the filtration associated with the continuous tensor product structure of the Fock space. The 
entries W\(t), . . . ,W m (t) of the vector W(t) are self-adjoint operators on which are associated with 
the annihilation and creation operator processes of external boson fields. Omitting the time argument, 
suppose the quantum Ito table of W is given by 

dWdW T := (dWjdW k )^ jMm = Qdt, (2) 

where Q. := ((Ojk)i^j,k^m is a constant complex positive semi-definite Hermitian matrix of order m. Its 
entrywise real part 

V =ReQ = (Q + Q T )/2 (3) 

is a positive semi-definite symmetric matrix since Q. = Q.* )p 0, with (•)* := ((-)) T the complex conjugate 
transpose. The imaginary part of the quantum Ito matrix Q. in © specifies the CCRs between the fields 
as 

[dW,dW T ] := {[dWj,dW k })^ ]Mm = Udt, J := 2ImQ, (4) 

where J G A m in view of the Hermiticity of Q.. The quantum Wiener process W represents a quantum 
noise which drives a QSDE 

dX = Fdt + GdW, (5) 
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governing the Heisenberg dynamics of the vector X of system variables. The rc-dimensional drift vector 
F(t) and the dispersion (n x m)-matrix G(t) of this QSDE are expressed as 

F:=i[H,X]+Sf(X), G:=-i[X,h T ] (6) 

in terms of a system Hamiltonian H, the quantum Ito matrix Q. of the process W and system-field 
coupling operators h\,...,h m which are assembled into a vector h := (h k )i^ k ^ m . Here, jSf denotes the 
Lindblad superoperator [15] whose action on the vector X is given by 

i m 

: = o L G>jk(hj[X,h k ] + [hj,X]h k ). (7) 

j,k=l 

Regardless of a particular form of H and h (which usually are functions of system variables), the CCR 
matrix in (OQ) is preserved in time due to the Heisenberg unitary evolution 

£(0 = ^(0^(0)^(0 (8) 

of observables £ on the system-field composite Hilbert space Jf? ® & (in particular, the lifted system 
variables X k (0) <g> /jr). Here, £/ (?) is a unitary operator on eg) with initial condition U (0) = /jf 
and (•)' is the operator adjoint. The quantum expectation of the observable © is defined as 

EC(0:=Tr(p(0)C(0), (9) 

where the initial system-bath density operator p (0) is assumed to be the tensor product 

p(0) := G7(0)®1) (10) 

of the initial plant state 07(0) on J^f and the vacuum state u := |0)(0| of the external fields associated 
with the vacuum vector |0) in where the Dirac bra-ket notation [20] is used. In the next section, 
we will specify the energetics of the system being considered. 

III. Hamiltonians with Lur'e type perturbations 

For what follows, the system-field coupling operators hi, . . . ,h m in © are assumed to be linear with 
respect to the system variables: 

h: = MX, (11) 

where M E W nxn is a constant matrix. That is, hj is a linear combination of Xi, . . . ,X n whose coefficients 
form the jth row of M. Then, in view of the CCRs (OQ) and the bilinearity of the commutator, the 
dispersion matrix G is a constant real matrix: 

G= ~i[X,X T ]M T = 0M T =.B. (12) 

Also, suppose the system Hamiltonian H in © is described by 

H:=Ho+j^(p k (Y k ), Y k :=c T k X, (13) 

k=l 

where 

H :=X T RX/2 (14) 

is a nominal Hamiltonian, specified by a real symmetric matrix R of order n (the space of such matrices 
is denoted by §„). Here, <pi, .. . ,<p s : R — >• K are continuously differentiable functions which are applied 
to self-adjoint operators Yi,...,Y s , assembled into a vector 

Y:=(Y k )i^ s = C T X, C:=[ci ... c J , (15) 



where the columns of the (n x s) -matrix C are the vectors c\,...,c s e R" from ([T3l . In combination 
with the linear system-field coupling (fTTT ). the quadratic Hamiltonian //o would lead to a linear quantum 
system, dynamically equivalent to the open quantum harmonic oscillator [4], [6] which is a basic model 
in linear quantum control [13], [22], [25]. Therefore, the functions <pi,...,<p 5 , which, in general, are 
not quadratic, can be interpreted as an unmodelled part of the Hamiltonian (fl3l) playing the role of a 
perturbation to the nominal quadratic Hamiltonian Hq. The following lemma specializes the QSDE © 
for the system under consideration. 

Lemma 1: For the open quantum system with the CCRs (OQ), linear system-field coupling operators 
(fTTI) and Hamiltonian (TTSl . the vector X of system variables satisfies the QSDE 

dX =Fdt + BdW. (16) 

Here, 

F=AX + ®CZ, A:=®R + BJM/2, (17) 
where the matrices J and B are defined by © and (fT2l . and 

Z:=(Z k ) i<k<s , Z k :=(p' k {Y k ), (18) 

with the derivatives <p' k of the functions <p k being evaluated at the operators Y k . 

Proof: Being a particular form of ©, the QSDE ([T6l is obtained by substituting (fTTT)— (IT3T) into 
©-(17]). The linear part AX of the drift vector F in ([171) comes from the quadratic Hamiltonian //o in 
(PT41) and the linear system-field coupling operators (fTTI ) as 

i[H ,X) = ®RX, Jzf (X) = BJMX/2, (19) 

whose derivation is well-known in one form or another [4]. The term ®CZ in (fTTI) originates from the 
relationship i[(p(Y k ),X] = i<p'(Y k ) \Y k ,X] = ®c k Z k , which follows from the commutator identities [20, pp. 
38-39], the notation and the CCR [Y k ,X] = -[X,Y k ] = -[X,X T ]c k = -i®c k between Y k and X in 
view of ©• ■ 
Therefore, the effect of quadratic perturbations described by 

<Pk(y) = Ykf/2, (20) 

with constant coefficients yi e R, is equivalent to modifying the nominal quadratic Hamiltonian 
Hq in (TTSl as R i-)- R + Cdiag 1 ^ yt ^ i ,(y-)C T , in which case, the drift vector F in (fT71) inherits the linear 
dependence on the system variables from ([T9l . If the functions <pi, . . . , <p s are not quadratic, F becomes 
nonlinear with respect to X. Such perturbations are similar to those in the classical Lur'e systems [16]. 
In Sections [VT - IVIIIL we will investigate the influence of the Lur'e type perturbations on the quadratic- 
exponential moments of the system variables. This study will be based on the more general results of 
the next section on exponential moments of adapted quantum processes. 

IV. Exponential moments of quantum processes 
Let be an adapted quantum process on the system-field product space Jf®^ satisfying a QSDE 

dt=fdt + g T dW, (21) 

driven by the quantum Wiener process W with the Ito matrix Q. from ©, where the drift / and the 
dispersion vector g := (gj)i^j^ m are adapted quantum processes. We assume that f(t), gi(t) , ■ ■ ■ , gm(t) 
are self-adjoint operators on if ®^ for any t ^ 0, and so is £(t). Such a QSDE can be obtained from 
©, or its specialization (fT6l-(fT8l, if ^ is a function (for example, a polynomial with real coefficients) 
of the system variables. In this case, bounds for the moments E(£ r ) of computed for positive integers 
r over the vacuum state of the external fields in the sense of ©, (flOl) . would guarantee a statistically 
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"localized" behaviour, that is, robust stability, of the quantum system [12], [26]. However, more subtle 
information on the system dynamics is provided by the exponential moment 

E(f) := Ee^ W = Tr(p(0)e^ W ). (22) 

Here, £ may be multiplied by a real, which, for simplicity, is not done in the present study since 
the required effect can be achieved by an appropriate scaling of the processes / and g in (l2TI) . The 
following Lemma |2] employs the ideas of [17], [18], [28, Appendix B], [31] on parameter differentiation 
of exponential operators to develop a representation for the quantum Ito differential of the exponential 
=Lr>o ( D r / r ^ see a l so [14, PP- 480-482], [27, pp. 200-206] on the exponentials of unbounded 
self-adjoint operators. A straightforward computation of is complicated by the noncommutativity 
between £ and d^. To formulate the lemma, let S x denote a linear superoperator, associated with <jj(f) 
(at an arbitrary moment of time t ^ 0) and acting on linear operators r\ on Jfl? ® & as 

g x (ri) :=e- A «ije*« = £ L^adjft) = e"^ ft). (23) 

Here, A is a real parameter, and ad^(rj) := [£,Tf] = — dx<§x{r\) \ X=Q is the negative of the infinitesimal 
generator of the superoperators with the group property SxoS^ = and identity element <£o- A 
particular case of the group property, S x <S-\ji = implies that 

e^(T]) = ^_ 1/2 (T])e^ 2 , (24) 

e^ftWO = 4-i/ 2 (T?)^-V2(0e l/2 (25) 

The definition (1231) also shows that the superoperator carries out a similarity transformation and is, 
therefore, homomorphic: 

S x {nQ = S x {i])S x {Q. (26) 



Lemma 2: The quantum Ito differential of the exponential of the adapted quantum process t, in (1271) 
is representable as 

de^ = S /2 (adt + p T dW)^ /2 . (27) 

Here, 

«:= / 1/2 fe(/)+^U) T ^/ A ^(g)d/i)<tt, (28) 

.7-1/2 V 7-1/2 / 

4/2 



j8 := (j8 7 ) KKm := / ^(g)dA (29) 

7-1/2 

are adapted quantum processes, self-adjoint at any time, and the superoperator S x from (|23l) applies to 
the vector g entrywise. 

Proof: For any positive integer r, repeated application of the quantum Ito rule d(r)Q = (dT])£ + 
7]d£ + (dr/)d£ to the rth power of £ yields 

+ £ ^(d£)$*(d£)^>-* 
where the second sum represents the Ito correction term. Hence, 

de ^ = y mil i y ^(djljMlj! . (30 ) 

jlk %o u+k+i+2)\ ■ 
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We will now use the following identities for Euler's multivariate Beta functions 

f (i-xy'i k M. = j\k\/u+k+i)\, (3d 

Jo 

I (l-A) J '(A- i u)VdAd i u = jlkl£\/(j + k + £ + 2)l, (32) 

J A 

which hold for all nonnegative integers j, k, £, with the integration in (l32l) being carried out over a 
planar simplex 

A:={(A,^)GR 2 : l^^^O}. (33) 



In view of (131b . the first sum in (f301 takes the form 



L TTTTjTiTT = L L ^ dA 



= / 1 e ( 1 - A ^(d^)e A ^dA=e^ / 1 c T A (d^)dA, (34) 
jo Jo 



where (1231) is used. The right-hand side of (1341 is recognizable as the Gateaux derivative of the 
exponential e^ [17, Eq. (10)], [28, Eqs. (B5), (B6)], [31, Eqs. (2.1), (4.1)], which, in our context, 
is evaluated in the direction of the Ito differential d^. Similarly, substitution of (l32l into the second 
sum in (l30l) yields 



(j + k + £ + 2)\ 

JU w\ dAdM 

e( 1 - A )^(d^)e^-^^(d^)e^dAd J u 



A 

J # A (d£K M (d£)dAdju 

-1 / M 



jf *xM)(f S^)d^dX, (35) 



where the multiple integral employs the structure of the simplex (1331) . This corresponds to the second 
derivative of the exponential [31, Eq. (11.6)]. By substituting (l34l) and (l35l) into ((301) . it follows that 



de^ =e^ jf 1 <f A (d£) (l +^ ^(d|)d/x)dA 

1 /2 A 

=e^ 2 / 4(d^)(l+/ # M (d£)d/iWe^ 2 . (36) 

7-1/2 V '-1/2 / 



'1/2 / /-A 

1/2 V ./-1/2 

Here, the second equality is obtained by using the identities (l24l) . (l25l) and appropriately translating the 
limits of integration. Since the Ito differential dW commutes with the adapted processes taken at the 
same or earlier moments of time, then 

S x (dS ) = $x (f)dt + S x (g) T dW, (37) 

where Sxis) a ^ so commutes with dW. By combining (1371 with the quantum product rules dtdW = 
and ©, it follows that 

^(d£K M (d£) = S k {gf€lS^g)dt. (38) 
6 



The representation (|27b is now obtained by substituting (|37l ), (|38l ) into (136*1) and using ([281) , (l29l ). Finally, 
the self-adjointness of the operators ce(?) and j8i (?),..., j8 m (?) can be verified directly from (128V (l29l 
by using the symmetry of the integration limits ±1/2 about the origin and the identity 

fab)? = *- X tf) (39) 

which follows from (ED and = £. Indeed, ([39]) implies that (^i(/)) t = and 
(£'h(g) T Q,£'n(g)y = £-p i (g) i: Q.£'_x(g), since the drift / and the entries of the dispersion vector g in 
((211) are self-adjoint operators and Q. = Q.*. Hence, 

. rl/2 ,f rl/2 rl/2 

/ * x {f)M) = / *_ A (/)<* = / * x (f)dX, 

W-l/2 / 7-1/2 7-1/2 

and 

( jT < r A (g) T fl<f M ( g )dAd/i) t = Jj?-^g) T &g_ x {g)&x&ii 

= f ^(g) T n^(g)dAd)U, 

which proves that both parts of a in (T2~8T) are self-adjoint, and a similar argument applies to the entries 
of /3 in d29l . Here, the set A := {(A,/l) G M 2 : 1/2 ^ A ^ /l ^ -1/2}, associated with (|33l . is invariant 
under the area-preserving linear transformation (A,/l) >->■ (— jU,— A). ■ 

If £ were a classical diffusion process, then d<^ would commute with thus making ad^(d(§) vanish 
and implying that (d^ ) = d^ . In this commutative case, the representation (1361 ) would reduce to 
de 1 ' = e 1 ' (di§ + (d^) 2 /2), leading to the Doleans-Dade exponential [3], with the 1/2-factor coming from 
the area of the simplex (|33l , and (d<^) 2 = g T Q.gdt being the quadratic variation of Therefore, the 
relation (|27T ) can be regarded as a noncommutative quantum counterpart to the classical stochastic 
exponential. The following theorem applies Lemma [2] to the dynamics of the exponential moment (l22l) . 

Theorem 1: The exponential moment (1221) of the quantum process t, from (|2TT) satisfies a differential 
equation 

S = E(e^W /2 ), (40) 

where the process a is defined by (T2~8T) . 

Proof: The commutativity between dW and adapted processes, combined with the product structure 
of the system-field density operator (flOl ). imply that 

E(e^ 2 /3 T dWe^ 2 ) = E(e^ 2 /3e^ 2 ) T EdW = 0. 

Therefore, since the noise term on the right-hand side of (|27l) does not contribute to S, the latter reduces 
to the average of the drift term e^ 2 ae^ 2 , thus proving (|40l . ■ 
In the context of the exponential moment dynamics, the quantum process a, which linearly enters 
the right-hand side of (|40l and is computed according to (|28l ). will be referred to as the rate process. 

V. Risk-sensitive moments of system variables 
As an adapted quantum process <§, we will now take a quadratic form of the system variables 

£:=X T lLY/2, (41) 

specified by a matrix IT G S n , so that £;(*) is a self-adjoint operator on Similarly to the usual 

quadratic forms of commuting variables, the diagonalization of IT allows t, to be represented as a linear 
combination of squared observables 
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whose coefficients <j\ ^ . . . ^ o n are the eigenvalues of II, with Vi , . . . , V n the corresponding orthonormal 
eigenvectors in W 1 . Therefore, if II )p 0, that is, 0\ ^ 0, the operator <§(f) is also positive semi-definite. 
The exponential moment (122b . associated with (l4Tb . takes the form 

S := Ee xTra / 2 (42) 

and will be referred to as the risk-sensitive moment (RSM) of the system variables. This extends the risk- 
sensitive storage functions from classical stochastic systems [19] to the quantum setting. In addition to 
the time dependence, the RSM S depends on II, which plays the role of a matrix-valued risk-sensitivity 
parameter. Its asymptotic behaviour for small II is described by 

S= l + (n,ReE(XX T ))/2 + o(n), 11^0, (43) 

where (K,L) := Tr(K*L) denotes the Frobenius inner product of real or complex matrices, and E(XX T ) = 
ReE(XX T ) + i®/2 is the matrix of second moments of the system variables. Since e 1 ' + % (this 
property is inherited from the exponential function on the real line due to the spectral theorem for 
self-adjoint operators), then S ^ E(/ + £) = 1 + (II,ReE(XX T ))/2, so that the first two terms on the 
right-hand side of (|43l provide a lower bound for S. By a slightly refined reasoning, if II )p 0, then 
e 1 ' )p E, r /r\, which yields an upper bound E(^ r ) ^ r!S for the higher-order moments of t, for all positive 
integers r. Now, to study the dynamics of the RSM S by the methods of Section [TV] we will first describe 
the evolution of the process <§. 

Lemma 3: For the open quantum system, specified by Lemma \T\ the process t, , defined by (14Tb . 
satisfies the QSDE (T2T)) whose drift / and dispersion vector g are computed as 

f=r + (X T (A T Tl + UA)X 

+x T n®cz-z T c T ®nx)/2, (44) 

t: = (BVB t ,U)/2, (45) 
g=B T UX, (46) 

where the matrix V is given by ©. 

Proof: By applying the quantum Ito formula to (|4T|) and using ©, (fT6b . it follows that 

2d^ = X T TldX + (dX) T UX + (dX) T TldX 
=X T Tl(Fdt + BdW) + (Fdt +BdW) T UX + dW 7 B T UBdW 

= (X T TIF+F T UX + (BVB T ,n))dt + 2X T TlBdW, (47) 

where Tr(^5 T IIB) = (BVB J ,11) by the antisymmetry of the matrix J in ©. It now remains to substitute 
(fTTT) into (gTJ) to verify that ^ is governed by (ED with ([44l)-(l46l ■ 
Application of Theorem [T] to the RSM (|42l) requires computation of the rate process oc, defined by 
(|28T) . for the drift / and the dispersion vector g of the QSDE (1271) which are specified by Lemma [3] 

VI. Computation of the rate process 

The process g in (14"6l) is linear with respect to X, whereas / in (144b is a quadratic function of the 
system variables perturbed by the terms which are bilinear in X and the Lur'e type nonlinearities Z from 
(fT8~b . Since the processes / and g enter the rate process a in (l28b as arguments of the superoperators S\, 
we need the following lemma which allows to be computed for linear, quadratic and more general 
nonlinear functions of the system variables by taking advantage of the CCRs (OQ). 

Lemma 4: The action of the superoperator S^, associated with the system variables by (|23b and (14Tb . 
on the vector X is equivalent to the left multiplication by a matrix exponential: 

£ X {X) = K X X, K x ;=j X&n = Kl x -, (48) 
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where is the CCR matrix of the system variables from (Q3. 

Proof: The superoperator ad,e applies to vectors entrywise and commutes with the left multi- 
plication by a complex matrix, so that ad^(/VX) = Nad^(X) for any N G C nxn . Hence, by using a 
straightforward induction, initialized with 

ad|(Z) = -j©riX (49) 

in view of (OQ) and (14TT) . it follows that ad^(X) = (— i&U) r X for any nonnegative integer r. Substitution 

of the last relation into yields S X (X) = £ l>0 (a©II) r X/r! = e' A0n X, which establishes ©. ■ 
Lemma |4] shows that the superoperator <£^, associated with (f4T|) . performs a Lie algebraic similarity 
transformation [31, Section 6] of the system variables. In particular, this implies that 

^(x l f(u T X)) = x i r(u T K x X) (50) 

for any vector u G C" and at least for analytic functions i/a of a complex variable. It suffices to verify 
(|5Q|) for monomials y(z) '■ = z r of arbitrary nonnegative integer degrees r. To this end, note that 

^((u T xy) = (^(u T x)Y = (u T ^(x)y = (u T K,xy\ 

which is obtained by combining the homomorphic property (|26l) with Lemma 0] Another identity, which 
follows from (|49l , is given by 

ad^ (XX T ) = ad^ (X)X T +Xad^ (X) T 

= i(xx T n®-®nxx T ), (51) 

where use is also made of the Leibniz product rule I] £] = [<!;, T]] £ + T] [£, £] f° r tne commutator. The 
fact that the right-hand side of (IFTT ) is quadratic with respect to X (as is £ in (|4T1)). is closely related 
to the property that the commutator of quadratic forms of annihilation and creation operators is also a 
quadratic form of these operators [26, Lemma 4]; see also [10, Appendix B] for other relevant results. 
Now, the linearity of g in (l46l) with respect to the system variables allows Lemma 0] to be applied as 

^ (g) = B T Yl<? x (X) = B T TIK X X. (52) 

For what follows, we assume that the system dimension n is even, and 

det0^O, n^O. (53) 

Then, in view of d^K^ = i®HKi, the representation (l52l implies that 

f X ^(g)d^i = B T n[ X KpdvX 

J -1/2 J -1/2 

=B r n(i®n)-\K x -K_ 1/2 )x 

= iM(K x -K_ l/2 )X, (54) 

where the matrix M = —B J ® 1 specifies the system- field coupling CCD). It will be convenient to extend 
the notation (K,L) to the case of a real or complex matrix K := (Kjk) and an appropriately dimensioned 
matrix L := (Lj k ) of operators on a common space by (K,L) :— Y.j,kKjkLjk- Substitution of (l52l) and 
(|54|) (which both are linear with respect to X) into the Ito correction part of the rate process a in (|28T) 
leads to a quadratic form of the system variables 

t 2 fe(g) T n t MsW)dX=X T TX = (T,XX 7 ), (55) 

7-1/2 V 7-1/2 / 
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where T is a complex Hermitian matrix of order n, which is computed as 

r := i f' 2 KjllBSlM(K k -K_ l/2 )dX. (56) 

In view of (|53l) . the matrix 011 = n _1 / 2 v / n©V / rTv / n is obtained by a similarity transformation from 
v^nQv^l G A„, with y/Tl a real positive definite symmetric matrix square root of IT. Hence, 011 is a 
diagonalizable matrix whose eigenvalues are pure imaginary and symmetric about the origin [8], that 
is, representable as =b'<0fc, where (Oi, . . . , (0 n / 2 are all real. In this case, the matrix exponential K x in (|48l ) 

is isospectral to a complex positive definite Hermitian matrix e ^^®^ with eigenvalues e ±Xc °k > 
(recall that the parameter X is real). Thus, 

K x =^Dx^\ (57) 
where the columns of the matrix *F e C nxn are the eigenvectors of 011, and 

Dx := diag (dj(X)) = diag (e ±A ^) (58) 

l</'<n Kifc<n/2 

is a diagonal matrix formed by the eigenvalues of Kx- We will now consider a linear operator Jf, which 
acts on an (n x n) -matrix P as 



pl/2 

X(P):= \ KxPK]dX=V(SO( x i> l PV T ))¥ T , (59) 

7-1/2 



where use is made of (iSTb . Here, denotes the Hadamard product of matrices, and *P T := (*P *) T . 
Also, the entries of the matrix 5 := (s jk)i<^j.k^n 6 §>„ are computed in terms of (158T ) as 

j ll \j{X)d k {X)dX. (60) 



and are all strictly positive. Therefore, the operator J£T in (1591) is invertible, and its inverse has a similar 
representation 

JT^/ 5 ) = x ¥(S e Q( x ¥" 1 P x l''' T )) x ¥ T , (61) 

where S e := (l/s/jfc)l<f' denotes the entrywise inverse of the matrix S. Also, the subspace of complex 
Hermitian matrices of order n is invariant under J£T. The significance of the operator is clarified by 
the relation 

■1/2 r \/2 



j <?x(XX T )dX=[ S x {X)Sx{XfdX 

7-1/2 7-1/2 
/•1/2 

/ i^XX^JdA = Jf(XX J ), (62) 

7-1/2 



1/2 

where use is made of the homomorphic property (|26l and Lemma 51 In the same vein, we define an 
(s x n) -matrix of operators 

•1/2 _ rl/2 

^(ZX T )dA= / , 

7-1/2 7-1/2 

where the rightmost equality employs the corollary d50l of Lemma @] applied to ([TBI as 

4(Z) = (^(cJ^X)) ia<s . (64) 



Q := I' 12 <?x(ZX J )dX = [ 1/2 ^x(Z)X T K T xdX, (63) 

7-1/2 7-1/2 
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Now, by substituting (|44l) . (1551) into (12~8l) , and using (1621 and (1631) , it follows that the rate process for 
the RSM <@2]) takes the form 

a =t + (A T n + riA, jt(xx t ))/2 + (r,xx T ) + (n, r)/2 
=t + (A T n + riA + utf-^T) , JT(XX T ) ) /2 
+ <n,r)/2. (65) 

Here, T is an (n x n) -matrix of operators associated with (|63) by 

T :=0C2-2 f C T 0. (66) 

Also, ^ := denotes the adjoint of the inverse operator J^~ l from (|61"| ) in the Hilbert space 

C' ?x ' 1 with the Frobenius inner product of matrices, which is computed as 

X^{P) = y-*(S e (^W))^ 1 . (67) 
VII. Dynamics of risk-sensitive moments 

Substitution of the rate process (|65l ) into (l40l of Theorem [T] yields the following equation for the 
RSM of the system 

d t S = tS + (A T n + ITA + 2 JT " f (f ) ,N))/2 

+ <n,E(e^ 2 re^ 2 ))/2. (68) 
Here, /V is a complex Hermitian matrix of order n defined by 

=JT(E(e^ /2 XxV /2 )). (69) 

The following lemma directly relates the matrix N with the RSM. 

Lemma 5: The matrix /V, defined by (1691) , is representable in terms of the RSM S from (1421 by a 
linear differential operator 

/V = 2dnS + j20/2, (70) 

where is the CCR matrix of the system variables from (H). 

Proof: Upon splitting XX T into the symmetric and antisymmetric parts as 

XX T = (XX T +(XX T ) T )/2 + i®/2, (71) 

the imaginary part of the matrix N in (|69l can be expressed in terms of (|421 as 

1 /2 

Im/V = E(e^ 2 y ^ ^(0/2)dAe^ /2 ) =20/2, (72) 

where use is made of the identity (®) = ® (recall that constants remain unchanged under the action 
of <%). Now, the process £ in (|4T|) depends linearly on the matrix IT G §„, and 

<9n^ = (^ T + (XX T ) T )/4 

coincides with the symmetric part of XX 1 in (|7TI) up to a factor of 2. Hence, a reasoning, similar to 
the proofs of Lemma |2] and Theorem [T| (except that no Ito correction terms arise in this case), allows 
the real part of /V in (l69l to be related to the Frechet derivative of the RSM S with respect to IT as 

1 /2 

Re/V = E(e^ 2 J ^^ X ((XX T + (XX 7 ) T )/2)dXe^ 2 ^ 

1 /2 

= 2E(e^ 2 j f i/2 ^(^)dAe^ 2 ) =2^S. (73) 
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The representation (1701 ) now follows from (1721 and (|73T ). ■ 
Lemma [5] allows the differential equation (1681 ) to be written in a more closed form. 
Theorem 2: Under the conditions (|53l) . the RSM (l42l) of the quantum system, described by Lemma [U 

satisfies a partial differential equation (PDE) 

d t Z =(A T Il + HA + 2T, dnZ) + (z+ (U, ®)/2)S 

+ (n,E(e^ 2 7e^ 2 )>/2. (74) 

Here, 

T:=ReJr _t (r), U :=ImJf- + (F) (75) 

are computed using (1561) . (l67l) . and T and T are defined by (|45l) . (l66l) . 

Theorem |2l which is the main result of the paper, shows that the unperturbed part of the RSM 
dynamics (the first line of (|74|) ) corresponds to a linear PDE of first order which can, in principle, be 
solved by the method of characteristics [5] using the solutions of the ordinary differential equations 
(ODEs) 

n+A T n+nA+2T = o. (76) 

This requires computation of the matrix-valued functions T, 13 of the risk-sensitivity parameter n y 
in (1751) and will be discussed elsewhere. The influence of the Lur'e type perturbations in the system 
Hamiltonian is described by the second line of (1741 ) and will be considered in the next section. 

VIII. Perturbations bounded using superpositive ordering 
Let L := (Ljk)i^.j,k^r be a matrix whose entries are linear operators on the system- field Hilbert space 
J$? '(g) J^", satisfying iJ := (L^i^^r = L, where, in application to matrices of operators, (-)^ := ((-)*) T 
denotes the transpose of the entrywise adjoint (•)*. Then for any r-dimensional complex vector u : = 
(uj) i^j^ r , the operator u*Lu :— Yfj k=\ u^jLjkUk is self-adjoint. We say that such a matrix L is superpositive 
if u*Lu y for any u E C. Since the superpositiveness extends the standard positive semi-definiteness 
(from a single operator to a matrix of operators), it will be written using the same symbol as L y 0. 
The superpositiveness of L is inherited by the matrix r\Lrf for any linear operator r\ on M*(£)&. Also, 
L y implies that EL y holds regardless of the reference density operator over which the expectation 
is computed. Similarly to the usual positive semi-definiteness, the superpositiveness induces a partial 
ordering among matrices with operator- valued entries. A straightforward example of a superpositive 
matrix is XX T . We will now use the concept of superpositiveness to specify a class of Lur'e type 
perturbations as those which satisfy 

T^oJ(r(XX T ), (77) 

where o > is a scalar parameter, and use is made of (l62l . (l66l . The matrix T in (f66l) depends on the 
risk-sensitivity parameter IT y through the matrix Q in (1631) since the matrix exponential K% in (|48T) 
involves II. The effect of this dependence on the description (1771 ) of the class of perturbations requires 
additional analysis which will be carried out elsewhere. To this end, we will only mention here that 
in a particular case of quadratic perturbations <pi,...,<p s given by (l20l) . the matrix Q becomes linearly 
related to the matrix on the right-hand of (1771 ) as 

Q= dmg(y k )C T Jr(XX T ). (78) 

Substitution of (1781) into (l66l) shows that the corresponding class (1771) of (quadratic) perturbations of 
the Hamiltonian can be described in terms of a Lyapunov operator P !->■ LP + PL T , where the matrix 
L := 0Cdiag 1 ^ yt ^ i (y^)C T does not depend on II. We will now discuss the role of the inequality (1771 ), 
which employs the concept of superpositiveness, in the dynamics of the RSM (l42l) for general (not 
necessarily quadratic) Lur'e type perturbations. 
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Theorem 3: Under the conditions of Theorem |2] and for the Lur'e type perturbations described by 
(1771) . the RSM of the quantum system satisfies a partial differential inequality (PDI) 

d t Z ^ (A T n + ITA + 2Y + an, <9nS) 

+ (t+(O,0)/2)S. (79) 
Proof: From (1771 ) and the properties of the superpositiveness, it follows that 

E(e^ /2 re^ 2 ) ^ c7E(e^ /2 jr(XX T )e^ /2 ) 

= oN = o(2d n E + iZ®/2), (80) 

where use is made of Lemma |5] Since n >- 0, then (|80l) leads to an upper bound for the perturbation 
term on the right-hand side of (I74T ): 

(n,E(e^ /2 re^ /2 ))/2 < a(Yl,2dnZ + iZ®/2} /2 

= a(n,dnZ), (81) 

where (11,0) = since II, belong to the orthogonal subspaces §„, A„. The PDI (|79l) is now obtained 
by combining fz!]) with <[8l]). ■ 
The PDI (|79l) can be treated as a multivariate partial differential case of the Gronwall-Bellman lemma. 
In fact, it is reduced to the standard univariate version of the lemma by noting that d t Z — (A T n + ITA + 
2Y + all, o\\Z) coincides with the total time derivative dZ/dt of the RSM E (as a function of t and n) 
along the trajectories of another characteristic ODE 

ri+A T n + IIA + 2Y + an = 0, (82) 

which is similar to (|76T ). Indeed, along the characteristics (|82l) . the PDI (|79l) is equivalent to an ordinary 
differential inequality 

dZ/dt^ (t+(O,0)/2)S, 
thus making the univariate Gronwall-Bellman lemma applicable to this case. 
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